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We consider the expected value for the total curvature of a random closed polygon. Numerical 
experiments have suggested that as the number of edges becomes large, the difference between the 
expected total curvature of a random closed polygon and a random open polygon with the same 
number of turning angles approaches a positive constant. We show that this is true for a natural class 
of probability measures on polygons, and give a formula for the constant in terms of the moments of 
the edgelength distribution. 

We then consider the symmetric measure on closed polygons of fixed total length constructed by 
Cantarella, Deguchi, and Shonkwiler. For this measure, we are able to prove that the expected value 
of total curvature for a closed n-gon is exactly f n + f 2 n- 3 • ^ s a consequence, we show that at 
least !/3 of fixed-length hexagons and i/n of fixed-length heptagons in R 3 are unknotted. 



1. INTRODUCTION 



The study of random polygons is a fascinating topic in geometric probability and statistical 
physics. Random polygons provide an effective model for long-chain polymers in solution under 
"^-conditions". There is an essential distinction between open polygonal "arms", which are easy to 
analyze because each edge is sampled independently, and closed polygons, where the closure con- 
straint imposes subtle global correlations between edges. The fixed-length open polygons with n 
edges in M 3 form a manifold which has a codimension 3 submanifold of closed polygons. To study 
random polygons, we must first fix a probability measure on open polygons and a corresponding 
codimension 3 Hausdorff measure on closed polygons. Given these measures, we can then study 
the statistical properties of the geometry and topology of polygons in each space. 

In recent work, two of us (Cantarella and Shonkwiler) ||2l presented a new measure on the 
space of closed ra-gons of fixed length constructed using a map from the Stiefel manifold V2(C n ) 
of orthonormal 2-frames in complex n-space to the space Poi3(n) of closed n-gons of length 2. We 
called this measure the symmetric measure. We computed exact expectations of radius of gyration 
and squared chord lengths with respect to the symmetric measure. These are global geometric 
invariants of polygons in the sense that they involve edges which are far apart along the polygon. 
In this paper, we are interested in the total curvature n, which is the sum of the turning angles at 
each vertex of the polygon. This is a local invariant of polygons, in the sense that it is determined 
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by pairs of adjacent edges. Though we will only consider total curvature in the present work, 
our methods should also apply to other local geometric invariants such as total torsion, which is 
determined by triples of adjacent edges. 

If we sample the edges in an open polygonal arm independently according to a spherically 
symmetric distribution, it is easy to see that the expected turning angle at a vertex is |. Thus the 
expected total curvature of an re-edge arm (which has re — 1 vertices) is |(n — 1). In a closed 
polygon, the edges are not independently sampled due to the closure constraint. Since the closure 
constraint involves all the edges, it is reasonable to expect that its effect at any given vertex should 
become negligible as n — > oo. Hence we expect that in a closed polygon, the expected turning 
angle at a vertex should also approach | . This is true, as we will see below. 

However, it is not true that lim^oo E(k, Pol 3 (n)) - |n = 0. In 2007, Plunkett et al. El 
numerically sampled random closed equilateral polygons of n edges and found that their total 
curvatures are equal to |n + a{n), where a(n) tended towards a constant near 1.2 as re — > oo. In 
2008, one of us (Grosberg) [4] presented an argument to explain why a(n) — > ^ = 1.1781 . . . 
for equilateral polygons. Numerical experiments O suggested that for the symmetric measure 
a{n) — > j. This raised an interesting question: why is the asymptotic value of a{n) different for 
the two measures? 

Theorem fT] answers this question. With respect to a probability measure /i satisfying mild 
hypotheses, the expected total curvature of a closed n-gon is |ra + a{n), where a(n) — > ^-^J- 
Here, mi and m,2 are first and second moments of edgelength. For random polygons sampled 
according to the symmetric measure of 0, this theorem shows that a(n) — > j, in agreement with 
our previous numerical experiment. 

The main result of this paper is that we are able to go much further for the symmetric measure 
and obtain an exact formula for the expectation of total curvature on Poi3(n). To do so, we show 
in Theorem [2] that the expectation of any scale-invariant function on P0I3 (re) in the symmetric 
measure is equal to the expectation of that function in a new measure called the Hopf-Gaussian 
measure. 

The Hopf-Gaussian measure is constructed by applying the Hopf map to a multivariate Gaussian 
distribution on quaternionic n-space. Since the coordinatewise Hopf map is a quadratic form on 
HP (cf. Section |3T| ), the coordinates of edges of polygons sampled according to this measure are 
differences of chi-squared variables. Hence, they have Bessel distributions. Using this fact, we 



determine the pdf of the sum of k edges in a random arm in Proposition 17 This enables us to find 



in Proposition 22 an explicit pdf for a pair of edges sampled from a random closed polygon. We 



then compute the expected value for turning angle by integration in Proposition 24 This gives us 



our main result (Theorem 15 1: the expectation of total curvature k for a random closed n-gon in 
the symmetric measure a is 

7r 7r 2re 
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This calculation gives us some new insight into these polygon spaces. For example, consider 
the old question: what fraction of the space of closed n-gons consists of knotted polygons? It has 
been proved that the fraction of unknots decreases exponentially quickly to zero in various models 
of random polygons |[3j HJ [15]]. There are also decades of computational experimentation on this 
question (cf. lfl2ll for references) which show that for small n, unknots are very common. Few 
theorems are known for specific values of n. Our total curvature theorem allows us to prove that, 
as measured by the symmetric measure, at least 1 /3 of the space of fixed-length hexagons in M 3 and 
i/li of the space of fixed-length heptagons in M 3 consists of unknotted polygons. 

These methods open up a number of new avenues for exploration and experimentation. The 
ability to write an explicit pdf for pairs, triplets, or other collections of edges raises the hope 
of computing expectations for other interesting scale-invariant functions, such as total torsion or 
average crossing number. Numerical integration with respect to these pdfs is also an effective 
method for approximating expected values. This can give significantly better results than averaging 
over large ensembles of polygons (cf. Section [5]). This will aid future research on these polygon 
spaces by allowing conjectured expectations to be tested to high accuracy. 

2. ASYMPTOTIC EXPECTED TOTAL CURVATURE OF POLYGONS 

The purpose of this section is to compute the asymptotic expected total curvature of closed 
random polygons under some reasonable hypotheses on the probability measure chosen for polygon 
space. We will denote the space of ra-edge open polygons (up to translation) in M. d by s£&(n) and 
the subspace of n-edge closed polygons in M. d by <%(n). Here we do not fix the lengths of the 
polygons, so «e^(n) consists of vectors of edges (el, ... , e^) G M. d x • • • x M. d = M. dn and ^(n) 
is the linear codimension d subspace of Arnid(ra) determined by the closure constraint ^ el = 0. 

We will say that a probability measure \i on £^(n) is generated by a spherically symmetric pdf 
g on W 1 when fi is the product measure g x • • • x g. The corresponding probability measure on 
S?dkn) is the subspace measure with respect to /x. Equivalently, we say that n is generated by g if 
arms are generated by sampling edges independently from g and closed polygons are generated by 
the same algorithm conditioned on closure. 

When the pdf g is spherically symmetric, we can write 

for some non-negative function /(r) so that f °° f(r) dr = 1 (and hence f Rd g{r) dVol^ = 1). 
In this case, the radial moments of g are the ordinary moments of /, and both are equal to the 
moments of edgelength with respect to /z. We denote these by 



m p := £(|e;M = E(\A P \9) = E(r p ;f). 



(2) 
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A number of standard probability measures on ^d(rt) are generated in this way. For instance, 
if f(r) = 5{r — 1), the resulting measure is the standard measure on ra-edge equilateral arms and 
the corresponding measure on ^(n) is the standard measure on closed equilateral polygons. If 
f(r) = Af(0, 1), the resulting measure on .sz^(n) is the standard measure on Gaussian random arms 
and the corresponding measure on 0P&{n) is the standard measure on Gaussian random polygons. 

If /i is generated by g, it is clear that the expected angle between two edges of a polygon in 
s^d{n) sampled according to \i is §, since the edges e~i are independently sampled from a spherically 
symmetric pdf on M. d . Thus, the expected value of total curvature on ^(n) is given by f (n — 1). 
Of course, an ra-edge closed polygon in ^(n) has an extra turning angle, so we might guess that 
the expectation of total curvature is |n instead. In fact, there is a curvature "surplus" in a closed 
polygon. We will now modify the argument in Hi to prove 

Theorem 1. For d > 2, if fx is a measure on s^dip) generated by a spherically symmetric pdf g 
which is bounded on IR d and has finite radial moments mi, vri2, and to 3 as in ([2]) and we take the 
corresponding subspace measure on &d( n )> then the expected value of total curvature on ^(n) 

approaches | n + gTP^i^Wiyi) ^ as n ~ * °°> wnere B is the Euler beta function. 
In particular, when d = 2 and 3 we have 

E{n;^ 2 {n),p )^-n+ x - and E(k; ^{n ~ -n + 1 - 

2 7T m-2 2 8 TO2 

for large n. 

Proof of Theorem^ We will assume for the duration of the proof that ^(n) (for any n) has a 
fixed probability measure p generated by a fixed spherically symmetric pdf g on M. d given by 
dif) = f(\?\)/ Vol S^^f^" 1 as in ([T]), and that ^(n) has the subspace measure induced by p. 

Let the Green's function Gfc(r) be the probability density of the end-to-end vector fin 
with respect to dVoV. We can write this explicitly as 

Gk(r) = J g{ei) ■ ■ ■ g(e k ) 8{ei H h e k - r ) dVol^ • • • dVol ?fc . 

If we consider the joint probability distribution of all edges in a closed polygon, we can treat it 
as a conditional probability on a set of edge vectors e\ , . . . , e n conditioned on the closure constraint 
^ = 0. This conditional probability can then be written as 

P(ei, . . . , e n ) dVol Sl ■ ■ ■ dVolg n = dVol^ • • • dVol^, 

where C n = G n (0) is the codimension d Hausdorff measure of the closed polygon space ^(n). 
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Given this joint distribution on all the edges, we can integrate out all but two of the edges to get 
the joint probability distribution on two consecutive edges, P(e*j, e$+i). Since this is independent 
of i, we may as well consider the case % = 1: 

P(ei, e 2 ) dVol el dVol e - 2 = g(ei)g(e 2 ) G "- 2( ~ Cl ~ 62) dVol* dVol e - 2 . (3) 

In other words, in order for the edges e\ and e 2 to come from a closed polygon, the remaining n — 2 
edges must connect the head of e 2 to the tail of §i. 

Finding exact expressions for G n - 2 and C n is quite challenging in general, but we can approx- 
imate both fairly easily by observing that the failure-to-close vector for an element of £fd(k) is 
just the sum of the edges. Given that we are sampling edges of our arms independently, that the 
third moment of g is finite, and that g is a bounded density on M d , the vector local limit theorem of 
Bikjalis CQ implies that the pdf of the (normalized) failure-to-close distribution converges in sup 
norm to the pdf of a normal distribution. 

We can recover the parameters of this normal distribution by noting that spherical symmetry 
implies that the mean of the failure-to-close distribution is zero, the variance of each coordinate of 
an edge vector is m 2/d, and the coordinates of an edge vector are uncorrelated. Therefore, the pdf 
of ^Gk{r) converges in sup norm to the pdf of the d-dimensional normal distribution 

N (0, diag (mafa . . . , ma/d)) , 

where diag(ai, . . . , etfc) is the diagonal k x k matrix with entries oi, . . . , a^. 

In particular, as n — > oo, we have that G n - 2 and C n = G n (0) are asymptotic in sup norm to 

/ d Y /2 ( -df 
G n - 2 (r) ~ — — exp 



2-7r(n — 2)777.2 / \2(n — 2)7712 
d 



d/ 2 



2irnm,2, 

where r = \r\. From Q, then, we see that the pdf P(e[, e$) is sup norm close to the function 

P(ei, e 2 ) ~ g(ei) 5 (e 2 ) ( exp ( ^±^1 ) . (4) 



77-27 V 2 ( n - 2 ) m 2 

Let 9(ei, e^) be the angle between e\ and e 2 , which is to say the turning angle between the two 



edges. We will now prove that E{6) | + B(i ^%l%% /2) ^ as n 



(d-l) B((d-l)/2,(d+l)/2) m 2 



00. 



First, for any e > we may choose p so that the integral of 6{e\, e*2)P(ei, e^) over the comple- 
ment of the ball B(p) of radius p centered at the origin obeys 

J 8(e u e 2 )P(e u e 2 ) dVol el dVol e - 2 < e. (5) 

W 2d -B(p) 
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To see this, observe that P(e±, e 2 ) is a pdf on M? d , so its improper integral over the entire space 
converges. This means that the L 1 norm of P on the complement of a ball of radius p goes to as 
p — > 00. But 9 is bounded by it so ||0(ei, e 2 )P(e*i, e* 2 ) 1 1 1 < 7r||-P(ei, e*2) 111 where the norms are 
over the complement of the ball B(p). Choosing p large enough that the rhs is less than e yields ([5]). 

Similarly, g{e[)g{e 2 ) is the pdf of a two-edge arm, so its improper integral over M 2d converges 
as well. Since g has finite first and second moments on M. d , this product has finite mixed moments 
of order up to 2 on M. 2d . In particular, for any quadratic polynomial q(ri,r 2 ) with coefficients 
bounded by ±A 2 we may choose p(A 2 ) so that we have 

J q{ri,r 2 ) g(ei)g(e 2 ) dVol^ dVol^ 2 < e. (6) 

R 2d -B(p) 



We now turn to the interior of the ball. Since the pair (el , 62) is in the interior of B{p) in R , we 

may choose n large enough that j^^m * s as c l° se to zero as we m particular, we may choose 
n large enough that the exponential in (|4|) is sup norm close to its linear Taylor approximation on the 

entire ball. Further, we can approximate ( J by 1 + ^ and ^2 by \ and we have P(e\, e^) 



sup norm close to the following function over the entire ball: 

(V) 



P(ei,e 2 ) ~ g{ei)g(e 2 
The expected value of Q{e\, e%) is just 
0(ex, e 2 )P(e*i, e 2 ) dVol^ dVol^ 



d d\ei\ 2 d\e 2 \ 2 d(ei,e 2 ) 



n 2nrri2 2nm 2 nm 2 



= J 9(ei , e 2 )P(el , e 2 ) dVol ?1 dVol e - 2 + J 6{ex , e 2 )P(ei , e 2 ) dVol ?1 dVol e - 2 . (8) 

B{p) M? d -B{p) 

By (|5]>, the second integral on the right is small and we can ignore it. Since the ball is a bounded 
domain, the fact that P(e*i , e 2 ) is sup norm close to the approximation in (|7]) tells us that the integral 
of the bounded function 9 against the approximation is close to the first integral on the right. 

Consider the approximation (|7]) to P(e*i , e 2 ). For large enough n, the quantity in square brackets 
is a quadratic polynomial in r\ and r 2 with coefficients between —2 and 2, so the inequality ((6]) 
applies with p > p(V2). Hence, since 9 is bounded, its integral against the approximation over 
the ball is close to its integral against the approximation over all of M 2rf . In other words, we can 
approximate the first integral on the rhs of ([8]> by the integral 

d d\e±\ 2 d\e 2 \ 2 d(ei,e 2 ) 



0(e 1 ,e 2 )g(ei)g(e 2 ) 



1 + 



n 2nm 2 2nm 2 nm 2 



dVole- dVol e - 2 
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We now evaluate the above integral. We will write e± and eb in spherical coordinates. Since 
the integrand is spherically symmetric, we can integrate out the angular coordinates of, say, e\ 
and assume that e\ lies along the z-axis. This produces a factor of VolS"* -1 . Since the rotation 
of e*2 hi the (d — l)-plane perpendicular to the e*i-axis does not change 9 or the approximation 
to P(e*i, e%), we can integrate out another VolS d ~ 2 . Since 9 is now the polar angle for e 2 and 



Vol S d ~ 
VolS d ~ 



r("/2) 



V / ^r((d-l)/ 2 ) 



, the integral reduces to 



r m 

nT ((d-l)/2) 



7T OO OO 




f(n) f(r 2 ) 

„<i-l d-1 







d dr\ dr\ 
n 2nm<i 2nm 2 



dr\r 2 
nm 2 



cos ( 



rf^rjf"" 1 sm d - 2 9 dn dr 2 d9. 



Since J °° /(n) drj = 1 and f °° rff(ri) drj = m p , integrating out r\ and r 2 yields 

nrri2 



((d-i)/ 2 ) J 
after integrating by parts. 



2 n(d - 1) B ((d-i)/2, (d+i)/2) m 2 



Since this is the expected value of the turning angle between two edges of the polygon, multi- 
plying by n yields the desired expression for expected total curvature. □ 



3. THE SYMMETRIC AND HOPF-GAUSSIAN MEASURES ON POLYGON SPACES 



Theorem [T] applies to a broad class of measures on polygon space, but not to certain highly 
symmetric measures defined in [2]. These symmetric measures on fixed-length polygons in space 
and in the plane are interesting for a number of reasons: they come from a natural geometric con- 
struction, expectations and moments of chordlengths and radii of gyration are exactly computable 
and scale like the corresponding expectations for equilateral polygons, and there is an algorithm 
for direct sampling from these measures which is fast (linear in the number of edges) and easy to 
code. Our goal for the rest of the paper is to determine the expected total curvature of polygons 
with respect to these measures. 

The symmetric measure is most naturally defined on the space of (open or closed) n-gons of 
fixed total length 2 in either R 3 or R 2 , which we denote by Arnid(n) for open polygons and Pol^ (n) 
for closed polygons in M. d . Of course we can extend the definition to the space of polygons of any 
fixed length by scaling. Viewed as a subspace of ^(ra), the space Pol^(n) differs from <5^(n) in 
that elements satisfy constraints on both closure and total length. Therefore, Theorem [T] does not 
apply to the symmetric measure. 
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However, since total curvature is a scale-invariant quantity and since ^(n) is a cone over 
Pold(n), the expected total curvature of polygons in ^(n) - which we can determine asymptot- 
ically using Theorem [T]- will be the same as the expected total curvature of polygons in Pold(n) 
provided that this expectation is computed with respect to a measure on ^(n) which is a prod- 
uct of some measure on the cone parameter and the symmetric measure on PoLi(n). Indeed, in 



Section 3.2 we will define the Hopf-Gaussian measure H on =e^(n) and ^(n) for d = 2, 3 by 



applying the Hopf map to the standard multivariate Gaussian measure on H™. The Hopf-Gaussian 

(ld-i n x <?, where x\d 



measure on stf d {n) turns out to be the product H = Xod-i x °"> where Xod-i is the chi-squared 



distribution with 2 d ~ 1 n degrees of freedom on the interval [0, +oo) which parametrizes the cone 
direction and a is the symmetric measure on Arm^n). Likewise, the Hopf-Gaussian measure on 
^(n) is the product H = x%i-\ x <T > where now a is the symmetric measure on PoLj(n). 

An immediate consequence of this construction is the following theorem, which is the central 
message of this section: 

Theorem 2. Suppose F : s^dip) — >• M is a scale-invariant function. Then the expected value of F 
over £/d( n ) with respect to the Hopf-Gaussian measure H is the same as the expected value of F 
over Armrf(n) with respect to the symmetric measure a; that is 

E{F;^ d (n),R)=E(F;Arm d {n),a). 

Likewise, if F : &d(n) —> R is scale-invariant, then 

E(F;0> d (n),R) = E(F;Pol d (n),a). 



As we will see, the Hopf-Gaussian measure satisfies the hypotheses of Theorem[T] so the com- 
bination of Theorems [T] and [2] will allow us to determine the expected asymptotic total curvature on 
PoLj(n) with respect to a from the first and second moments of edgelength on =sa^(n) with respect 
to H, which we compute in Section 3.3 These asymptotic total curvature expectations are given 
by: 

Corollary 3. For d £ {2, 3} and large n, the expected total curvature on PoLj(n) with respect to 
the symmetric measure is 

IT 2 7T 7T 

E(k; Pol 2 (n), a) ~ -n + -, E(k; Pol 3 (n), a) ~ -n + -. 

2 7T Z 4 



The value of ? for the total curvature surplus of polygons in M 3 agrees with our numerical 
experiments in . 

Of course, Theorem [2] applies to any scale-invariant functional on polygons, not just total cur- 
vature. We expect that it will be useful for determining the expected values of other interesting 
quantities such as total torsion and average crossing number. 
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3.1. Quaternions and the Symmetric Measure on Polygon Spaces 

In this subsection we recall the construction of the symmetric measure on Arm t j(n) and Pol^n) 
from [2J. Recall that these are spaces of arms and polygons of fixed total length 2. In principle 
everything could be scaled to any desired fixed length, but the choice of length 2 will be the most 
convenient. Since the translation of the following definitions and results to any other scale is 
straightforward, we will not discuss this scaling further. 

Definition 4. For d E {2, 3}, let Armd(n) be the moduli space of n-edge polygonal arms (which 
may not be closed) of length 2 up to translation in M d . An element of Arnid(n) is a list of edge 
vectors ei, . . . , e n € Mr whose lengths sum to 2. 

Consider the Hopf map from the division algebra of quaternions H to the space of imaginary 
quaternions (which we identify with IR 3 ) given by 

Hopf(g) = qiq, 

where q is the quaternionic conjugate of q. In coordinates, if q = (go, qi, (fe; (ft)* then 

Hopf (g) = (<?o + qf ~ ql ~ ql, 2g x g 2 - 2q q 3 , 2g g 2 + 2qiq 3 ). (9) 

We extend the Hopf map coordinatewise to a map Hopf : BP — > M 3n . Then Hopf is a smooth map 
from the sphere 5 ,4n ~ 1 of radius v2 in H n onto Arni3(n). Specifically, for q G S An ~ l the edge set 
of the polygon Hopf (q) is 

(ei, . . . ,e„) := (Hopf(gi), . . . ,Hopf(g n )). 

We call S^' 1 the model space for Arni3(re). 

Similarly, the restriction of Hopf to the 1 © j and i © k planes gives a map to the i © k plane, 
which we identify with C = M 2 . Specifically, 

Hopf (a + 6j) = i(a + 6j) 2 
Hopf(ai + 6k) = i(a + 6j) 2 . 

In other words, if z = a + 6j, then Hopf (z) = iz 2 and Hopf (iz) = iz 2 . Extending this map 
coordinatewise yields a smooth, surjective map Hopf : S 2n ~ l U5 2n_1 — > Arm2(n); consequently, 
the disjoint union 5 2n_1 U S 2n ~ l is the model space for Arni2(n). 

Definition 5. For d G {2, 3}, let Pol t j(n) be the moduli space of closed n-gons of length 2 up to 
translation in M d . 
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Since PoLi(n) C Arm^n), the inverse image Hopf (PoLi(n)) is well-defined and will be the 
model space for Pol^(n). To describe this model space for d = 3 as a subset of S^ n , the sphere 
of radius v2 in BP, it is convenient to write S An ~ l as the join S 2 "^ 1 * S 211 ^ 1 . In coordinates, the 
join map is given by 



where u, v 6 C n lie on the unit sphere and 6 G [0, "/2]. The Stiefel manifold V2(C n ) of Hermitian 
orthonormal 2-frames (u, v) in C n can be identified with the subspace 



and, as Hausmann and Knutson first observed 0, this manifold is precisely the model space for 
Pol 3 (n). 

Proposition 6 ([5]). The coordinatewise Hopfmap takes V2(C n ) C C x C = i" onto Pols(n). 

The key to proving the above proposition is to note that the Hopf map applied to a quaternion q 
can be written more simply by letting q = a + b] for a, b G C: 

Hopf (q) = Hopf(a + 6j) = (|a| 2 - \b\ 2 , 23?(a&), 2M(ab)) = i(\a\ 2 - \b\ 2 + 2abj). 

Let V2(W l ) be the real Stiefel manifold of orthonormal 2-frames in R n , which sits naturally in 
M n © j]R n . A result analogous to the above holds for planar polygons: 

Proposition 7 ((21). The coordinatewise Hopfmap takes V2(M. n ) U iV^M") onto Pol2(n). 

With these maps in place, we can define probability measures on the arm and polygon spaces by 
pushing forward measures on the model spaces. Since the model spaces are homogeneous spaces, 
it is natural to push forward Haar measure on the model spaces, which is what we did in [2]. Since 
Haar measure is also the measure defined by the standard Riemannian metrics on these spaces, this 
gives the following definition of the symmetric measure a on P0I3 (n) : 



The symmetric measures on the other arm and polygon spaces are defined analogously. 

The space Pol^n) is topologically the union of the spaces of polygons with fixed edgelengths 
n, . . . ,r n such that r i = 2, so any expectation over Pol^(n) with respect to the symmetric 
measure is a weighted average of the expectations over these spaces. In future work we intend to 
determine how the average is weighted and with respect to which measure on the fixed edgelength 
spaces. For equilateral polygons the answer is simple and pleasant: the restriction of the symmetric 
measure to the subspace of equilateral polygons is just the natural measure on this space, namely 
the subspace measure on n-tuples of vectors in the round S 2 which sum to zero. 




{(u,v^/4):(u,v)=0}cS' 



4n-l 
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3.2. The Hopf-Gaussian Measure on Polygon Spaces 



We would like to compute the expectation of total curvature over Pol^(n) with respect to the 
symmetric measure. Unfortunately, since elements of PoLi(n) are chosen from ^(n) by con- 
ditioning on the polygon both being closed and having total length 2, we cannot directly apply 
Theorem[T] However, since the total curvature k is scale-invariant and since ^(n) is a cone over 
PoLj(n) (with the cone direction parametrized by the length of the polygon), we have that 



for any measure p on ^(n) such that p = p x a for some measure p on [0, +oo). 

At the level of model spaces, the picture is clearer. Hopf maps HP onto 3/3 (n) and the image 
of the sphere of radius r is exactly the copy of Amid(ra) consisting of polygonal arms with total 
length r 2 . Moreover, the measure on this scaled copy of Arm^(n) is exactly the pushforward of the 
standard measure on the sphere of radius r, so it is the symmetric measure defined in the previous 
section. Therefore, we can define a measure on s^{n) which is the product of some measure 
on [0, +00) and the symmetric measure on Arni3(ra) simply by pushing forward any spherically 
symmetric measure on HP. Of course, if we may choose any spherically symmetric measure, 
the obvious choice is the multivariate Gaussian measure: this is both a spherically symmetric 
measure and a product measure on the coordinates and we can expect that the fact that the individual 
coordinate distributions are Gaussian will simplify our computations considerably. 

In the case of planar polygons, the model space for s$2 (n) is the explicit copy of C n U C n given 
by (IP © jlP) U (ilP © kIP), so we will push forward the Gaussian measure on C n : 

Definition 8. If 7 4n is the standard Gaussian measure on HP = M 4n , then the Hopf-Gaussian 
measure H on =2/3 (n) is defined by 



Likewise, if 7 2 ™ is the measure on (l n © jlP) U (ilP © kK n ) naturally induced by the standard 
Gaussian measure on C n = M? n , then the Hopf-Gaussian measure H on 3/2 (n) is defined by 



E(k; & d {n) ,p) = E(k; PoLj(n) , a) 





The fact that the multivariate Gaussian is a product measure implies that: 
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Proposition 9. The Hopf-Gaussian measure on ^(n) is generated by the pdf 

e -W/2 
9{rl = 16^ 

<2«<f Hopf-Gaussian measure on ^(n) is generated by the pdf 

g-W/2 

Proof. Since 7 4ra is a product measure on EP, its restriction to each H 1 factor is the standard four- 
dimensional Gaussian. In particular, for q = (qi, . . . , q n ) 6 W 1 sampled according to 7 4n , the 
qi S H are independent, identically distributed, and spherically symmetric. Therefore, the edges of 
the polygon Hopf (q) = (Hopf(gi), . . . , Hopf (q n )) are independent, identically distributed, and, 
since the Hopf map is SU (2)-equivariant, spherically symmetric. Therefore, the Hopf-Gaussian 
measure H on .2/3 (n) is generated by a spherically symmetric distribution on M 3 , which we now 
determine. 

The four real components of each qi are themselves Gaussian-distributed. Therefore, since 
I Hopf((7j)| = I qi 1 2 , each edgelength of a Hopf-Gaussian polygon is given by the sum of the squares 
of four Gaussian real numbers, so these edgelengths follow a chi-squared distribution with four 
degrees of freedom. Thus, the pdf of the edgelength distribution is 

— r /'2 

fir) = — r . 

But then the pdf g of the spherically symmetric edge distribution is 

1 e - |r1 / 2 

9{f) = iw /(|r1) = th^t 

as desired. 

The proof in the 2-dimensional case is completely parallel. □ 

To identify the model space for ^(n), note that, as with the Pol^(n) and Armd(n) spaces, 
we have that ^(n) C ^(n). Therefore, Hopf _1 (^(n)) is a well-defined subset of H n which 
will be the model space for ^(n). We can identify this model space more explicitly as follows. 
Focusing on the case d = 3 for the moment, since W 1 is the cone over g 4 "- 1 it is convenient to 
write W l as the cone of the join 5 2n_1 * S 2n ~ l . In coordinates, the "cone-join" map is given by 

(u, v, 9, s) 1 — y s(cos9u + sin#tTj) 

where n, v € C n are unit vectors, 9 £ [0, w /2] and s G [0, +00). The cone CV2(C n ) over the Stiefel 
manifold V^C") can then be identified with the subspace 

{(u,v,*/4,s) : (u,v) = 0} C H n 

and the proof of Proposition [6] generalizes to show that CV2(C n ) is the model space for ^(n): 
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Proposition 10. The coordinatewise Hopfmap takes CV2(C n ) C W 1 onto &^(n). 

If CV2(M. n ) is the cone over the real Stiefel manifold V^M"), then the same reasoning yields 
the analogue of Proposition [7j 

Proposition 11. The coordinatewise Hopfmap takes CV 2 (M. n ) U iCV 2 (M n ) C BP onto @> 2 (ri). 

Since &d{ n ) C =2^(n), we can define the Hopf-Gaussian measure on ^(n) as the subspace 
measure inherited from the Hopf-Gaussian measure on =sz^(n) from Definition [8] 

To prove Theorem [2j which says that the expected value of any scale-invariant function on 
polygons is the same whether we compute it with respect to the Hopf-Gaussian measure or the 
symmetric measure, it suffices to show that the Hopf-Gaussian measure is a product measure: 

Proposition 12. Suppose d=2or3. Then the Hopf-Gaussian measure H on <e/d( n ) is the product 
X2d-i n x °~ of the chi-squared distribution with 2 d ~ 1 n degrees of freedom on the interval [0, +oo) 
and the symmetric measure on Arm^(n). 

Likewise, the Hopf-Gaussian measure on &d{ n ) is the product X^d-i n x °~ of the chi-squared 
distribution with 2 d ~ 1 n degrees of freedom on [0, +oo) and the symmetric measure on Pol^(n). 

Proof. Since the Gaussian measure on HP = M 4n is S'0(4n)-equivariant, its restriction to the 
sphere S 4n_1 (r) of radius r is, after normalization, just the uniform probability measure on the 
sphere. Since Hopf (S' 4n_1 (r)) is the space of arms of total length r 2 , this means that the restric- 
tion of the Hopf-Gaussian measure on =2/3 (n) to this space is, after normalization, the symmetric 
probability measure a defined in Section [3~Tj Likewise, the restriction of the Hopf-Gaussian mea- 
sure on ^(n) to planar arms of total length r 2 is just the symmetric measure. 

Therefore, the measure on £^(n) is the product p x a for some measure p on the interval 
[0, +00), so it suffices to see that p is the chi-squared distribution. Since the interval parametrizes 
the total length of a polygon, we need to analyze the distribution of total length of polygonal arms. 
For q = (qi, . . . , q n ) S HP, the arm Hopf (q) G =2/3 (n) has total length 

^2 iHopf(<7i)i = ^2 = ^2 i^i 2 - 

Since each \qi | 2 is the sum of the squares of four standard Gaussians, the total length of the polygon 
Hopf (q) follows the standard chi-squared distribution with An degrees of freedom. Therefore, the 
measure p on [0, +00) is the measure induced by the chi-squared distribution with 4n degrees of 
freedom. Likewise, for polygonal arms in the plane, the measure on total length is induced by the 
standard chi-squared distribution with 2n degrees of freedom, since in that case each \qi\ 2 is the 
sum of the squares of two standard Gaussians. 

The fact that the Hopf-Gaussian measure on ^(n) is the product of the chi-squared distribution 
on [0, +00) and the symmetric measure on PoLj(n) then follows immediately from the definition 
of the Hopf-Gaussian measure on ^(n) C =2^2 (n) as the subspace measure and the fact that the 
symmetric measure on Polrf(n) C Arm^(n) is the subspace measure. □ 
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Theorem [2] now follows since a scale-invariant function is by definition independent of the first 
factor in the product decomposition of the Hopf-Gaussian measure. 



3.3. Moments of Edgelength, Expected Chordlengths, and Expected Gyradius 

By Proposition [9] the Hopf-Gaussian measure on s^d{ n ) f° r d £ {2,3} is generated by the 
spherically symmetric pdf 

fi -M/a 

g(r) 



4 d - 1 7r|f| ' 

This is certainly a bounded density with finite first, second, and third moments, so we can use 
Theorem[T]to compute the asymptotic expected total curvature on ^(n) with respect to the Hopf- 
Gaussian measure. To do so, we just need to know the first and second moments of edgelength. 
In fact, as we saw in the proof of Proposition |5J the edgelength distribution on «e^(n) is the chi- 
squared distribution with 2 d ~ 1 degrees of freedom, so the moments of edgelength are just the 
well-known moments of this distribution: 

Proposition 13. The pth moment of edgelength on &?d( n ) i s given by 

E{\e^\ ^(n),H) = Fpl, £(|e^; #/ 3 (n), H) = 2"(p + 1)! 

Note that the expected values of edgelength are 

E(\ei\;^ 2 {n),H) = 2, £(|e-|; *f 3 (n),H) = 4 (10) 
and the expected squared edgelengths are 

and E(\ei\ 2 ;£/ 2 {n),B) = 8, E(\ei\ 2 ;£/ 3 {n),R) = 24. (11) 

Using the above values for mi and 771,2, Theorem [T] implies that the asymptotic expected total 
curvature on ^(n) is given by 

£(K;^(n),H) ^ £n + - %^ 3 (n),H) ~ T^+T- 

Z 7T Z 4 



Since total curvature is scale-invariant, Theorem [2] implies Corollary |3j which says that for large 

7T 2 7T 7T 

E(k; Pol 2 (n), a) ~ -77 + -, E(k; Pol 3 (n), a) ~ -77 + -. 

Z 7T Z 4 



77 



Also, we can now compute the expected value of chordlength and radius of gyration for =2^(77) 
using results from 0. 
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Corollary 14. The expected value of the squared length of a chord skipping k edges on <K^(n) is 

J B(Chord(/c); s 6/ 2 (ra),H) = 8k, £(Chord(A;); s/ 3 (n), H) = 24fc. 
The expected squared radius of gyration for arms in £/d( n ) is 

E(Gyradius; */ 2 (n), H) = - K £(Gyradius; ^ 3 (n), H) = 4- 1 r J 



3n+l n+1 

Proof. Since our measure on «e^(rt) is invariant under rearrangement of edges, we can easily com- 
pute expected squared chord length and radius of gyration using Propositions 5.3, 6.3, and 6.5 
from |2]. Those propositions imply that 

£(Chord(A0; £? d (n), H) = kE(\ei\ 2 ; £? d (n), H) (12) 



and 

6(n+l) 

Substituting the second moment of edgelength from ( fTT| ) into (12i and ( [T3] > yields the desired 
results. □ 



£(Gyradius;^(n),H) = f ^" j ^ ) ^(|e;| 2 ; ^(n), H). (13) 



4. EXPECTED TOTAL CURVATURE OF POLYGONS WITH THE SYMMETRIC MEASURE 

Corollary [3] gave the asymptotic expected total curvatures with respect to the symmetric mea- 
sures as 

7T 2 7T 7T 

E(k\ Pol 2 (n), a) ~ -n + -, E(k; Pol 3 (n), a) ~ -n + -. 

2 7T 2 4 

Our aim in this section is to use the special properties of the Hopf-Gaussian measure to carry out the 
argument of Theorem[T]for space polygons with no approximations. This will yield the following 
exact expectation which holds for any n > 3: 

Theorem 15. The expected total curvature with respect to the symmetric measure on P0I3 (n) and 
the Hopf-Gaussian measure on (n) is given by 

E{K;Pol 3 (n),a) = E(k; ^ 3 (n),H) = -n + 



2 4 2n - 3 

Since total curvature is scale-invariant, the first equality is a consequence of Theorem[2j proving 
the second equality is the main task of this section. Also, since all triangles have total curvature 



2ir, Theorem 15 is trivially true for n = 3. Therefore, throughout the rest of the section we will 
assume n > 3. 

We will shortly be evaluating many definite integrals involving the Bessel function K u (z); to 
do so we will repeatedly avail ourselves of the following: 
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Lemma 16. For real \i > \v\ and a > 0, 



x 



e- ax KJax) dx = ^ ' ' K \ , '- (14) 

v ' 2^ r(/i + i/2) 



and 



7o « M M + ^ r(M/2 + l/2) 

Proof. Both equations follow from the identity (71 6.621(3)] 

7o X 6 dX " (a + r(/x+l/2) 2 'V ' 2'^ 2'a + /3y'' 

which holds for any complex numbers fj,,v,a,f3 with > |3 f ?(i / )| and 5ft(a + /3) > 0. The 

function 2F1 is Gauss's hypergeometric function. 

To get ( fl4| ), we simply use the fact that 2-^1(0, b; c; 0) = 1 for any a, 6, c; to get ( fl~5| ) we use 
Kummer's identity: 2 -Pi(a, 6; a - 6 + 1; -1) = rja+Slr^/^i) • D 



4.1. The Sum offc Edges in s^{n) 

Since the goal is to carry out the strategy from Section [2] with no approximations, we first need 
to explicitly determine the Green's function Gk- 

Proposition 17. The probability distribution of the vector f joining the ends of a k-edge sub-arm 
of an arm in i^(n) with the Hopf -Gaussian measure is spherically symmetric in M 3 and given by 
the following explicit formula: 

where r = \r\. 

Proof. Suppose q = (q\, . . . , q n ) 6 M k is sampled from the standard Gaussian distribution. Writ- 
ing qi = ai + bii + cy + d^k for each i = 1, . . . , k, the failure-to-close vector for the k-edge arm 
Hopf(g) € ^(n) is 

y^Hopf(gj) = Hopf ( a j + 6 j i + a j + d fc k) = ^^(af +&f -c 2 -df 1 2biC i -2a i di,2aiCi+2b i di) . 

Again, the fact that this vector follows a spherically symmetric distribution is a consequence of the 
fact that the Hopf map is S[/(2)-equivariant. 
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Since the aj, &j, Cj, d\ are chosen from standard (real) Gaussian distributions, the distribution of 
the projection of the failure-to-close vector onto the first coordinate is clearly the difference of two 
chi-squared distributions, each with 2k degrees of freedom. This is known to have a pdf given by 
a Bessel function distribution (8j Chapter 12, Section 4.4] in the form 

fM = J$m Kk -' lAm - (16) 

It is worth noting that this pdf was proved by McLeish lPTTTl to be the pdf of a product of a gamma 
(k, 2) variable and an independent standard normal. McLeish also works out the moments and 
cumulants of the distribution. 



Next, we will use a result of Lord IflOl which relates the pdf of the projection of a spherically 
symmetric distribution to the pdf of the full distribution. 

Lemma 18 ( 11101 Eq. (29)]). Suppose p(r) is a spherically symmetric distribution on M 3 and that 
the projection of p{f) to any radial line through the origin has the pdfp\{r), where r = \r\. Then 
p{r) is given by 

p(r) = -— p\(r) dVolf. 
2irr 



Using the pdf of the projection given in ( fT6] ), Lemma [T8] implies that the failure-to-close distri- 
bution is 

This can be re- written using the derivative identity K' u (u) = —^(K u _i(u) + K v +\(u)) and the 

recurrence relation K v {u) = K v _ 2 {u) + ^=^-K^i(u) (cf. [6, 10.29.1] for both) to get the 
desired expression 

r k - 3 / 2 K k _ a/2 (r/ 2 ) 
G*«dVoV= ^ dVol„ 

where r = \r\. □ 

to the 



Using the identity K_y 2 (z) = z (cf. [6, 10.39.2]) and specializing Proposition 17 

case k = 1, we see that the distribution of an edge in a Hopf-Gaussian arm is 

— r ji 

GAr) dVol r -= dVol r -. 

167rr 

Note that this is, as it should be, the same distribution for edges that we found in Proposition [9] 
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4.2. The pdf of Edges in 5% (n) 



Proposition 19. 77ze codimension 3 Hausdorff measure of ^(n) in ^(n) f/ie va/«e G n (0), 
which is given by 

r(n-3/ 2 ) 



Proof. We saw in Proposition 17 that 

2n _ 2 n _ 3/2 ^ , , ^ 

G n (r) = 3 " 3/2 W j . (17) 

7T / 2 1 (n) 

We must be careful evaluating this formula at r = 0, since the Bessel function K n _ 3 / 2 has a pole 
at 0. To rewrite (17) in a form which allows us to easily evaluate at r = 0, we use the general 
formula for Bessel functions of half-integer order from ||T6l p. 80, formula (12)]: 



/ \ /7r\ 1/2 (n + i)l 



\2z) ^ i\(n - i)\{2zY 

Writing n — 3 /2 = (n — 2) + 1/2, we see that G n (r) simplifies to 



2- 2( " +1) e"-/ 2 ^ (i + (n-2))! r(w _ 2) . 



vrr(n) E »!((„_ 2) -i)! 



In this form, it is clear that the only term in the sum which is nonzero at r = is the i = n — 2 
term and 

r _ 2~ 2{ - n+1 \2n - 4)! _ 2- 2 (™ +1 )r(2n - 3) 
n(0) ~ 7iT(n)(n-2)! ~ 7rr(n)r(n - 1) ' 

We can simplify this a bit further using the duplication formula for gamma functions [6, 5.5.5] 

T{2z) = Ti- 1 l 2 2 2z - l T{z)T(z + 1/2) to get 

^-^-sSr (18) 

as desired. □ 



We can use the pdf for sums of edges in (n) to write down the pdfs for single edges and pairs 
of edges in «^(n). 
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Proposition 20. The pdf of a single edge e*j in (n) with respect to dVol^ is spherically sym- 
metric on M 3 and given by the following function of Yj = \ei\: 



n — 1 



- r i/2 n-T/2 



2 2n - 2 7iT(n - 3/ 2 ) 
Proof. The pdf of a single edge is just 



r'r "K n _ 5/2 (n/2) . 



p (ei) = yr9{ei)G n -i(-ei); 

i.e., the probability of the ith edge being e*i and the remaining n — 1 edges summing to —Si, 
conditioned on the assumption that all n edges sum to zero. Since g{ei) = G\(ei), we can use 
Proposition 17 and Proposition 19 to arrive at the stated expression. □ 

Corollary 21. The moments of edgelength for polygons in ^(n) are 

E(n-, aV„),H) = ""^f;"^" 3 ' B(P + 2, n - 2). 

ant/ hence the expectation of squared chordlength and radius of gyration are 

'n-k\ 12fc(2n-3) 



S(Chord(A:);^3(n),H) 
S(Gyradius;^ 3 (n),H) 



n 

n — 1 



n + 1 
(2n-3). 



Proof. Using the pdf from Proposition 20 the expected value of |e*j| p = r- is 



rfP(ei) dVol e - 



OO f-TT p2-K 



rfP^&i) rf sin (f> d9 d(f) drj 



'0 J jo 

Writing out P(e*j) and integrating with respect to 6 and gives the pth moment 

U " 1 ^ - ri/ Vr +p - 3/2 iv- n _ 5/2 (n/ 2 ) dr, 

n-1 V^ r (P + 2 ) r ( 2ra + P- 3 ) 



2 2 ™- 4 r(n - 3/2) y 



2 2m " 4 r(n - 3/2) 



_ (n- l)r(2n + p-3) 
~ 2T(2n - 4) 



B(p + 2,n-2) 



using Lemma 16 and the duplication rule for the gamma function. 
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Since the measure on ^(n) is invariant under rearrangement of edges, we can use Proposi- 
tions 6.3 and 6.5 from [2] as in Corollary [T4| Specifically, 

£(Chord(fc);^ 3 (n),H) = f^j) kE(\ei\ 2 ; <? 3 (n),R) 



and 
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£(Gyradius; ^(n),H) = ( ^ ) £(|el| z ; ^b(n),H) 



By the first part of the proposition we have that 

£(|e-| 2 ;^ 3 (n),H) 



2. ^ _ 12(n-l)(2n-3) 



(n + 1) 



so the given formulas for E'(Chord(A;); ^(n), H) and .E(Gyradius; ^3(n),H) are immediate. 

□ 

Proposition 22. The probability distribution of a pair of edges el and el in ^(n) is invariant 
under the diagonal action of SO (3) on the pair of edges and invariant under rotations which fix 
one edge. Hence, it depends only on the lengths r\ and r 2 of the two edges and the angle 9 between 
them. It is given by the formula 

P(r 1 ,r 2 ,0)dr 1 dr 2 d6 = -r l r 2 e^ {ri+r2 h n -lK ri _7 ( -) sin 9 dri dr 2 60, (19) 

4- v /7rr(2n — 4) 2 V2/ 

where z = | e\ + e 2 \ = \Jr\ + r\ + 2rir 2 cos 0. 

Proof. As we saw above, the general form for the probability distribution of a pair of edges in a 
closed polygon in ^3 (n) is 

P(ei,e 2 ) dVol^ dVol^ = — dVol^ dVolg- 2 , (20) 

that is, the probability of the first two edges being el and el and the remaining edges summing 
to —el — el conditioned on the assumption that all n edges sum to zero. We computed Gk(r) in 
Proposition 17 and C n in Proposition 19 Using the formula for g(e) from Proposition |5J we get 

P(ei,el) dVol e - a dVola = - ^ ^(n^n-?/^ ^ dvd 

327T ' 2 rir 2 l (2ra — 4) 



where again z = \ e% + e 2 1 = \Jr\ + r| + 2rir 2 cos We can rewrite el in spherical coordi- 
nates el = (ri,4>i,9\). For each el, we can fix el as the z-axis of spherical coordinates for 
el = (r 2 , $2, 9 2 ). Here 9 2 is equal to 9, the angle between el and el. Observing that we can inte- 
grate out 4>i and (f> 2 immediately to get a factor of An 2 and 9\ to get a factor of 2, and recording the 
volume form in these coordinates as r\r\ sin 9 dri dr 2 d9 gives us the formula in the statement of 
the proposition. □ 
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Corollary 23. The pairwise distribution of edges e\ and e 2 in ^3(11) for n > 3 may also be written 
more simply in terms of the variables x = (r\ + f2)/2, y = [r\ — r 2 )/2, and z = \e.\ + e*2 1- In 
these variables, the probability distribution is given by: 

P(x, y, z) dx dy dz = ^^M_ e -z z n~s/ 2 R ^ ^ dx &y ^ (21) 

Proof. Computing the Jacobian of the map (r±, r 2 , 6) \-t (x, y, z), we see that its inverse determi- 
nant is 

2 csc(0) y/r\ + r\ + 2r\ri cos 9 
nr 2 



We then multiply the function in Proposition 22 by this determinant and substitute to obtain the 
statement of the Corollary. □ 

We now explicitly check that the total integral of the pairwise pdf is equal to 1 . This computation 
will serve as a warm-up for the more difficult definite integrals ahead. 



Since y does not appear in the pdf P(x, y, z) in (21_), we will integrate with respect to y first. 
The triangle inequality states that r\ < r 2 + z, so y < z /2, and similarly r 2 < r\ + z, so — z /2 < y. 
We will next integrate by x, since x does not appear in a Bessel function. There is no upper bound 
on x = (ri + r 2 )/2, but since r\ + r 2 > z, we know x > z /2. We will integrate with respect to z 
last, and here the limits are simply and oo. Thus, we are trying to show that 



T(n) 



2^T{2n - 4) J 
The y and x integrals are simple, and leave us with 



[■CO pOO pz/2 

/ / / e- x z n - 5 / 2 K n _ 7/2 ( z h) dydxdz = l. (22) 

JO Jz/2 J-z/2 



w2- 4) r^""K^ mi , 



(23) 



by Lemma 16 This check gives us confidence that our pairwise pdf is correct so far. 



4.3. Expected Turning Angles and Total Curvature 



Proposition 24. The expected value of the turning angle 9 for a single pair of edges in ^3 (n) is 
given by the formula 
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Proof. Our overall strategy will be to use the formula for P(x, y, z) from Corollary [23] to write this 
expected value as 

poo poo P z /2 

= / / / 0(x,y, z)P(x,y, z) dy dx dz. 

JO Jz/2 J-z/2 

We begin by writing the turning angle 9 in terms of the x, y, and z variables as 

' z 2 - 2{x 2 + y 2 ) 



9(x, y, z) = arccos 



2(x 2 -y 2 ) 



Since P(x, y, z) does not depend on y, the first integral is accomplished by integrating this function 
with respect to y. Integrating by parts with dv = 1, we get 

t){x,y,z)dy = irz+ / = dy. 

-*/a J-*/* (y 2 - x 2 )\J z 2 - Ay 2 

Now, making the trig substitution y = | sin ^ yields 



/ r /2 sin 2 ^ / Z^ 2 / 8x 2 \ 

nz + V^ 2 ~ z 2 sinV _^ = vr, + y^ 2 - z 2 ^ [I - 8x2 _ z2 + z2cos27p ) W 



by using the identity sin 2 ip = 1 c ™ 2 ^ - This is simply 

ttz + its/ Ax 2 — z 2 — 2irx 
using the indefinite integral identity [7, 2.558(4)] 

dip 2 (a — b) tan f/2 

: arctan ■ 



a + 6 cos ip Va 2 - 6 2 Va 2 - 6 2 

Since the function P(», y, z) is in the form e~ x f(z), we must now do the pair of integrals 

fOO rOO 

tt / e~ x (z - 2x) dx + 7r / e~V^ 2 - z 2 dx. 



The first integral is simple and has the value — 2ire z / 2 . To do the second integral, we make the 
change of variables x = { z /2)t to get 



poo 2 poo 

/ e~V4x 2 - z 2 dx = — e~*3 (i 2 - 1)3 dt = ztfi , 

A/2 2 A 
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recognizing the last integral as a form of the integral representation for Bessel K v (6l 10.32.8] that 
we have used before. We have now shown that 

m = 2 ff- > 4) (^°°-2e-* / '3"-'" /2 K„_ V3 (»/2) dz + J" z"-'"K, K„_,/ a (f2) dz 



where the integral is, as usual, computed using Lemma 16 



The remaining integral in ( |25| ) is more interesting, as it involves a product of Bessel functions. 
Making the substitution z = 2u, we can then use the Nicholson integral representation for the 
product of Bessel functions [6, 10.32.17] to rewrite the integral as 



2r(2n - 4) 



z n ~^K x (*/2)K n _ 7/2 (*/a) dz 



Con-l/2p/ \ t-co / poo \ 

T(2n-4) I cosh (( 9 /2-^)^y o n"- 3 / 2 ^ n _ 5/2 (2 U cosht)dnJ dt. 

Again, the inner integral is a power of z multiplied by a Bessel function of az (here a = 2 cosh t) 
and can hence be evaluated using Lemma [16] We now have the integral 

^ 2 5 / 2 -™(n-2)(n-l) / cosh ((9/ 2 - n) t) sech™" 1 / 2 * dt. (26) 
■/ o 

We will integrate this using the general integration formula 

/ cosh ax secrr cc ax = W , (27) 

7o V 2 r(p) 

where Pjf is the associated Legendre function. This is valid when /3 — a > and /3 + a > 0. It is 
a specialization of [7, 3.5.17]. We set a = 9 /2 — n and /3 = n — 1/2 and see that /3 — a = 2n — 5 
(which is positive since we have assumed n > 3) and /3 + a = 4. Applying the formula shows that 
the integral of ( |26] ) is equal to 

37r 3/2 2 3-n (n _ 2 )(n - l)P 4 1 r n "(0)r(2n - 5) = 4(n-2)(n-l) 

r(n-i/2) (2n-5)(2n-3) 7F ' 

where we have used the general formula for the value at zero of associated Legendre functions 
given in (6j 14.5.1] and the duplication formula for gamma functions to simplify the form on the 
left hand side. Combining ( |25] ) and ( [28] ), we see that 

TI — 1 7T 7T 7T7T2 

^-aSTj'-a + terB-S + S^ (29) 

as desired. □ 
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Theorem 15 now follows by multiplying by n. We get an interesting corollary of this theorem 
for n = 6 and n = 7; since the expected value of total curvature is less than An, some hexagons 
and heptagons must have total curvature less than 47r and hence be unknotted by the Fary-Milnor 
theorem |pT3l . 

Corollary 25. If we measure volume using the symmetric measure on polygon space, at least 1 /z 
of the polygons in Pol3(6) and 1 /u of the polygons in Pol3(7) are unknotted. 

Proof. Let x be the fraction of polygons in Poi3(n) with total curvature greater than 47r. By the 
Fary-Milnor theorem, these are the only polygons which may be knotted. We know that any closed 
polygon has total curvature at least 2ir, so the expected value of total curvature satisfies 

E(k; Pol 3 (n), a) > Airx + 2vr(l - x). 

Solving for x and using Theorem [15] we see that 

(n-2)(n-3) 



x < 



2(2n - 3) 



For n > 7, this bound is not an improvement on the trivial bound x < 1, but for n = 6, we get 

x < 2 /3 and for n = 7, we get x < 10 /n, as desired. □ 

If instead we had let x be the fraction of polygons in Pols(n) with total curvature greater than 

2ttB, then x < (g^)2(2n-3) anc ^ tn ^ s §i ves constraints on the fraction of knots with bridge number 
B. For example, when n = 8, 9, 10, or 11 this gives us some information on the fraction of 
34)ridge knots since there are 34)ridge knots (e.g. 819) with stick number 8. 



5. NEW NUMERICAL METHODS FOR RANDOM POLYGONS 



In Q, we gave a fast sampling algorithm for random polygons in Poi3(n) which is guaranteed 
to sample directly from the symmetric probability measure on this space given a supply of nor- 
mal random variates. The ensembles of polygons generated by this algorithm are as good as the 
underlying ensembles of normals, so the quality of this sampling algorithm cannot be improved. 
However, when one is computing the expected value of a geometric functional on polygon space 
such as total curvature, averaging over a large ensemble of sample polygons is simply Monte Carlo 
integration over a very high-dimensional space. The numerical accuracy of this method is neces- 
sarily limited. 

The framework above yields a much better method for computing expected values of functions 
like total curvature which are both scale-invariant and sums of quantities that are defined locally 
on a given polygon: integrate directly against the pdf for a finite collection of edges in a closed 
n-gon in ^(n). Carrying out this method is not trivial if the integrand, like turning angle, has 



25 



a singularity in most coordinate systems. However, the results are worth it. Figure [T] shows the 
curvature "surplus" term of total curvature minus |n plotted with data from sampling and the 
number of correct digits obtained by sampling and by numerical integration. 
Total Curvature Surplus Correct Digits 

10 




Number of edges n 



15 20 
Number of edges n 

FIG. 1: On the left, we see the average total curvature surplus n — plotted for ensembles of 5 million 
polygons with 5 to 25 edges generated by the direct sampling algorithm of |2| (dots), plotted together with 

(curve) and the asymptotic value of this 
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the exact formula for this surplus of j 2 „" 3 gi ven by Theorem 
surplus of j given by Theorem [T] (line). The fact that these three computations are in agreement serves 
as a useful check on our work above. On the right, we see the number of correct digits in the calculation 
of average total curvature by careful numerical integration of turning angle against the pairwise pdf of 



Corollary 23 (top line, about 10 correct digits), averaging over ensembles of 5 million polygons (middle 
line, about 4 correct digits), and averaging over ensembles of 1 million polygons (bottom line, about 3 
correct digits). We can see that we obtain significantly better results by numerical integration. 



6. FUTURE DIRECTIONS 



It is clear that the methods above have many interesting applications. For instance, we can hope 
to compute the expected total curvature for plane polygons as well. In this case Theorem [T] shows 
that the expected curvature surplus is 2 /tt. However, the integral seems somewhat forbidding and 
we do not have an explicit conjecture for a closed form. 

More promising is the direction of extending our results to other functionals on space polygons. 
Since Theorem [2] tells us that the expected value of any scale invariant functional over the space of 
Hopf-Gaussian polygons is equal to the expected value over polygons with the symmetric measure, 
we can expect to compute a number of other interesting expectations this way. For example, the 
argument of Section [2] can certainly be generalized to predict an asymptotic expected total torsion 
of \n — j. It would be interesting to find the exact expectation of total torsion. It seems possible 
to apply these methods to average crossing number as well, which is certainly a topic for further 
investigation! 
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In Corollary [25] we used our expectation for total curvature with respect to the symmetric mea- 
sure to bound the fraction of unknotted fixed-length hexagons below by 1/3 and the fraction of 
unknotted fixed-length heptagons below by 1/11. Such bounds are clearly too small: numerically 
sampling ensembles of 5 million polygons shows the fraction of unknots to be roughly 9999/io,ooo 
for hexagons and 2499 /2500 for heptagons 1 . Our bounds even significantly underestimate the frac- 
tion of fixed-length hexagons and heptagons with total curvature less than Air, which a similar 
experiment with 5 million samples reveals to be approximately 91.4% and 63%, respectively. We 
could improve our bounds by computing the variance of total curvature, or even by finding an 
explicit expression for the total curvature pdf. However, this will involve a more subtle global 
analysis of the correlations between turning angles in closed polygons, and we leave this topic for 
future work. 
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